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Abstract. The closely related issues of (1) reconstructing density and 
velocity fields from redshift survey and distance indicator data and (2) 
determining cosmological parameters from such data sets are discussed. 
I focus on possible explanations for the discrepant values of j3 = Q 0,6 /6 
resulting from analyses that employ various reconstruction methods. Al- 
though no firm conclusions are reached, possible routes towards resolution 
of these discrepancies are suggested by the discussion. 



1. Introduction 

The original charge to our panel, which consisted of myself, Marc Davis, Avishai 
Dekel, and Ed Shaya, was to discuss "Reconstruction Methods." To narrow the 
discussion, we decided in advance to concentrate on the problem of deriving j3 
from comparison of redshift survey and distance indicator data. This seemed 
appropriate given the panel membership: all of us have worked on this problem, 
and together we span the spectrum of results, from low values of (5 suggestive 
of a low density {£Im ~ 0.2) universe to high values consistent with Qm = 1- 

While it may appear that we have strayed from the theme of "reconstruc- 
tions," /3-determination is, in fact, closely linked with how we reconstruct the 
underlying velocity and density fields from the observable data. I clarify this 
point in § 2 below. In § 3, I review recent results for (5 from peculiar velocities 
and consider possible explanations for their wide variance. I conclude in § 5 
with a brief look to the future of the subject. 



2. Reconstructions and [3 

Our theoretical formulation of cosmic dynamics involves the mass overdensity 
field S p (f) and the peculiar velocity field v p (f). Linear gravitational instability 
theory tells us that these fields obey the simple relation V • v p = —f(^M)^p, 
where f(Q M ) ~ V% 6 . 

Such continuous fields are abstract mathematical constructs; reality is more 
complex. Discrete entities — galaxies — sparsely populate the universe, and we 
measure their redshifts and estimate their distances. Reconstruction Methods 
are the numerical algorithms we employ to turn these real-world data into rep- 
resentations of the underlying fields S p and v p . 

Several obstacles stand between the measurements and successful recon- 
struction of the density and velocity fields. First, there is biasing: how is the 
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galaxy overdensity 5 g related to the desired <L? In the oft-used if oversimplistic 
linear biasing paradigm 5 g = b5 p , and the linear velocity-density relation be- 
comes V • v p = —f35 g . To the extent linear biasing and linear dynamics hold, 
only (3, not Q^i is measureable by this method. Incorporating "realistic" mod- 
els of biasing — and it is by no means clear that such models exist yet — presents 
another challenge to reconstruction methods. Reconstruction of v p is if anything 
an even more daunting task: we estimate only the radial component of this field, 
at discrete, irregularly distributed positions, with large errors that grow linearly 
with distance. 

We cannot invoke the theoretical velocity-density relation, and thus measure 
/3, until we have reconstructed the underlying fields. If our reconstructions of 
Sg{f) and/or v p (f) are flawed, we'll get the wrong value of (3. And given the 
difficulties such reconstructions face, errors in (3 are, perhaps, inevitable at this 
early stage in our understanding. 

Figure |l] summarizes the relationship between reconstruction methods and 
/^-determination. The box at the bottom concerns the nature of the comparison 
by which (3 is obtained. If one reconstructs v p from the peculiar velocity data 
and compares its divergence with 5 g from redshift survey data, one is doing a 
density- density (d-d) comparison. Alternatively one can reconstruct v p (f) from 
redshift survey data (for an assumed value of (3), using the integral form of 
the linear- velocity density relation, and compare with observed radial peculiar 
velocity estimates. This is known as the velocity-velocity (v-v) comparison. The 
distinctions between the two approaches are subtle but important, as discussed 
below. 

3. Discrepant (3 Values, and Possible Explanations for Them 

It was not until the early 1990s that measuring (3 by comparing peculiar ve- 
locity and redshift survey data became a realistic goal. The advent of full-sky 
redshift surveys — notably those based on the IRAS point source catalog — and 
large, homogeneous sets of Tully-Fisher (TF) and related data were the key de- 
velopments. The earliest attempt was that of Dekel et al. (1993), the so-called 
"POTIRAS" comparison. In this procedure the velocity field is reconstructed 
using the POTENT algorithm, and its divergence compared with the galaxy den- 
sity field from IRAS (a d-d comparison). Dekel et al. (1993) found /3i = 1.29,[] a 
value widely taken as support for the then-popular Einstein-de Sitter paradigm. 
The POTIRAS analysis has since been redone with much improved peculiar ve- 
locity data, with the result f3j = 0.89 ± 0.1 (Sigad et al. 1998), still consistent 
with a critical density universe unless IRAS galaxies are strongly anti-biased. 

Between the first and second POTIRAS papers, a number of studies, based 
on the v-v comparison, arrived at markedly lower values of (3j. Shaya, Tully, 
& Peebles (1995) estimated M = 0.17 ± 0.1 (/?/ « 0.35 if IRAS galaxies are 
unbiased) using the Least Action Principle to predict peculiar velocities. Willick 
et al. (1997) and Willick & Strauss (1998) used the VELMOD method to find 



x Here and below we apply the subscript I to f3 when the galaxy density field is obtained from 
IRAS. Other redshift surveys yield different overdensities and thus different /3s. 
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Figure 1. A schematic representation of the relationship between 
input data, reconstructions, and derivation of cosmological parameters. 



f3j = 0.5 ±0.06. The above results were obtained from TF data; an application of 
similar methods to more accurate (but far more sparse) supernova data yielded 
f3j = 0.4 ±0.1 (Riess et al. 1998). (This list is not exhaustive, and apologies are 
due to authors not cited; my aim is illustrative rather than comprehensive.) 

Thus, it appears that a somewhat bimodal distribution of (5 values has 
emerged. The d-d, POTIRAS method produces (3i close to unity; the v-v com- 
parisons, several based on the same redshift and velocity samples as POTIRAS, 
yield fli ~ 0.5. Neither the d-d nor the v-v comparison is inherently more valid; 
both are firmly grounded in linear gravitational instability theory. What, then, 
could be the cause of the discrepancies? 

Neither the panelists nor the conference participants arrived at a satistfac- 
tory answer to this question. However, the following list of possible explanations 
proved to be fertile ground for discussion. The list is accompanied by my own 
(biased) commentary on the salience of each explanation. 



• Malmquist bias. This famous statistical effect has long been invoked as a 
cause of discrepant conclusions in astronomy. In fact, it is now a non-issue. 
Methods are now available for rendering Malmquist and related biases 
inconsequential; see Strauss & Willick (1995) for a detailed discussion. 
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• Nonuniverality of the TF and other distance-indicator relations. 

Another non-issue. It is easy to make the general argument — nonuniversal 
distance indicators imply spurious peculiar velocities — but there is little 
or no evidence implicating such effects in the problem. 

• Calibration Errors in Peculiar Velocity Datasets. Such errors may 
indeed be present; new results from the Shellflow project (see the paper 
by Courteau in these proceedings) suggest that the widely used Mark III 
catalog has across-the-sky calibration errors. However, these errors affect 
mainly bulk flow estimates, not the value of (3i (see Willick &z Strauss 1998 
for a clear demonstration of this). 

• Non-trivial biasing. It is now widely believed that biasing is not only 
nonlinear, but stochastic and nonlocal (see the panel review by Strauss 
in these proceedings). If so, it is perhaps unsurprising that different ap- 
proaches produce different results, as they may in fact be measuring differ- 
ent things. This argument may well contain a kernel of truth, but I am not 
convinced it is the whole story. I was impressed by the work of Berlind, 
Narayanan, and Weinberg on this subject (see their contributions to the 
proceedings), which shows that for all but the most contrived models of 
biasing, the value of (3 one obtains is relatively insensitive to methodology. 

• The density-density versus the velocity-velocity comparison. This, 
I think, is the central issue. We may be asking too much of our distance 
indicator data when we use them to derive the full 3-D velocity field and its 
derivatives, as is required for the d-d comparison. In the v-v comparison, 
by contrast, the distance indicator data is used essentially in its raw form; 
only the redshift survey data, which is intrinsically more accurate, is sub- 
ject to complex, model-dependent manipulation. From this perspective, I 
would argue that the low values, (3i s=s 0.4-0.5, that have come out of the 
v-v analyses are more likely to be correct. 

4. A Quick Look to the Future 

If we learned anything from our panel discussion, it was the usual lesson: bet- 
ter data will help. Fortunately, some already exist, and more are on the way. 
The Surface Brightness Fluctuation (SBF) data of Tonry and collaborators (see 
Tonry's paper in these proceedings) are now available for comparison with the 
IRAS redshift data; preliminary results are reported by Blakeslee in these pro- 
ceedings. SBF distances are considerably more accurate than TF distances and 
promise a much higher-resolution look at the velocity field. Adam Riess reported 
new, and also very accurate, SN la distances for nearby galaxies. The extant TF 
data are not likely to increase dramatically in the short term, but will be recali- 
brated by the Shellflow program, and perhaps the different TF data sets (Mark 
III, SFI, Tully's Catalog) will be merged into a larger, homogeneous catalog. 
In the longer term (~ 5-10 years), they will be supplanted by much larger and 
more uniform TF data sets that will emerge from the wide-field infrared surveys 
currently under way, 2MASS and DENIS (see the contributions by Huchra and 
Mamon in these proceedings). 
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I draw yet another lesson from our panel discussion: we have not yet fully 
developed the analytical methods needed to deal with nonlinearities in the uni- 
verse, both dynamical and biasing-related. In this regard I would encourage 
theorists to address the problem of, how, given a particular biasing model, do 
we predict the peculiar velocity field from redshift survey data, taking full ac- 
count of nonlinear dynamics? When reliable methods are developed for doing 
this, and tested against N-body simulations, we will be able to fully exploit the 
improved peculiar velocity data sets of the future. 

Acknowledgments. I thank my fellow panel members, Marc Davis, Avishai 
Dekel, and Ed Shaya, for stimulating discussions. 
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